clear all
**Save to update later
global BS_fold_old = "$BS_fold"
global working_old = "$working"
global tempfold_old = "$tempfold"
global clean_old = "$clean"

global spec = "new_klm"
local spec = "$spec"

global tempfold = "$BS_fold/local10/" // **UPDATING TEMPFOLD TO BE A LOCAL FOLDER FOR EACH BOOTSTRAP;
global workingnew = "$BS_fold/local10/inputs/"
global BSfoldnew = "$BS_fold/local10/final/"
global cleannew = "$BS_fold/local10/inputs/"

use "$working/prowess_wits_estimation_qonly_klm.dta", clear
save "$workingnew/prowess_wits_estimation_qonly_klm.dta", replace

use "$working/first_last_years.dta", clear
save "$workingnew/first_last_years.dta", replace

use "$working/wholesale_priceindex.dta", clear
save "$workingnew/wholesale_priceindex.dta", replace

global BS_fold = "$BSfoldnew" //  **UPDATING BS_fold TO BE A LOCAL FOLDER FOR EACH BOOTSTRAP;
global working = "$workingnew" //  **UPDATING working TO BE A LOCAL FOLDER FOR EACH BOOTSTRAP;

global Zbase = "lag_l lag_m lag_D_ln_micro_inv lag_D_micro_inv lag_I_sq_micro_inv lag_I_sq_micro_inv_D lag_I_sq_lm lag_I_sq_micro_inv_km lag_I_sq_micro_inv_D_km"
global ZnonX_base = "lag_D_ln_micro_inv lag_D_micro_inv lag_I_sq_micro_inv lag_I_sq_micro_inv_D lag_I_sq_micro_inv_km lag_I_sq_micro_inv_D_km"
global INS =  "dm_lag_q_m dm_lag_l dm_lag_m dm_lag_D_ln_micro_inv dm_lag_D_micro_inv dm_lag_I_sq_micro_inv dm_lag_I_sq_micro_inv_D  dm_lag_I_sq_lm dm_lag_I_sq_micro_inv_km dm_lag_I_sq_micro_inv_D_km"
local spec = "new_klm"


********************POINT ESTIMATES********************************************
do "$code/estimation/flexible/prod_est_qonly_klm_opf.do"

global saveoutputs = "$working"
do "$code/estimation/flexible/D_est_TFPQIV_Qklm.do"

do "$code/estimation/flexible/data_work_Q_klm.do" 
global datafiles = "decomposition_SS_ProdQ_new_klm" 
do "$code/estimation/flexible/crosscheck.do"

global datafiles = "decomposition_SS_ProdQ_new_klm" 
global samplename = "ProdQ"
do "$code/estimation/flexible/reg_shell_short.do" 

********************POINT ESTIMATES********************************************



***GRAB BS SAMPLES
do "$code/estimation/grabBS/10_Q_klm_grabBS.do"

***RUN BS
local bsamplesize = 100

forvalues bs = 1(1)`bsamplesize' {
	*******Decomposition Sample*******
	global bs = `bs'
	
	clear 
	use "$BS_fold_old/record/Qklm/bssample_`bs'"
	
	save "$tempfold/bssample1", replace
	*Sample for second step (firm year)
	bysort co_code1 bs_co_code1 year: generate count_thing = _n
	keep if count_thing==1
	keep co_code1 bs_co_code1 year
	save "$tempfold/bssample2", replace

	do "$code/estimation/flexible/prod_est_qonly_klm_opf_bs.do" 
	
	global saveoutputs = "$tempfold"
	do "$code/estimation/flexible/D_est_TFPQIV_Qklm_bs.do"
	
	global saveoutputs = "$tempfold"
	do "$code/estimation/flexible/data_work_Q_klm.do"
	
	global datafiles = "decomposition_SS_ProdQ_new_klm" 
	do "$code/estimation/flexible/reg_shell_short_bs.do"  
	
}

global BS_fold = "$BS_fold_old" //  **UPDATING BS_fold TO BE THE MAIN PATH AGAIN;
global working = "$working_old" //  **UPDATING working TO BE THE MAIN PATH AGAIN;
global tempfold = "$tempfold_old"
global clean = "$clean_old"
